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A concentration difference of particles across a membráne perforated by pores will induce a diffu- 
sive flux. If the diffusing objects are of tlie same length scale as the tlie pores, diffusion may not be 
simple, objects can move into the pore in a configuration that requires them to back up in order to 
continue forward. A configuration that blocks fiow through the pore may be statistically preferred, 
an attracting metastable state of the systém. This effect is purely kinetic, and not dependent on 
potentials, friction or dissipation. We discuss several geometries which generate this effect, and 
introduce a heuristic model which captures the qualitative features. 



PACS numbers: 87.16.dp, 66.10.cg, 83.10.Rs 

Managing flows of ions and molecules across mem- 
branes is a centrál task of living systems. Celí mcmbrancs 
contain a wide variety of pores, which can sort, drive, 
and rectify fluxes of different particles in both directions. 
There are many models for the controlled leaks, and one- 
way flows of particular ionic species, which are useful to 
a celí [Q. The simplest leak is passive diffusion through 
a pore, driven by a concentration gradient. If we assume 
nonintcracting point particles, the flow can be described 
by Fick's law. The flux is linear in the concentration 
gradient across the membráne, the diffusion coefficient is 
independent of concentration. However, in a biological 
context, diffusing objects are of finite size, often closely 
packed, and often not small compared to a pore diameter. 
These conditions can produce diffusion with a highly non- 
linear dependence on concentration, and lead the systém 
to states far from the linear Onsager regime of nonequi- 
librium thermodynamics. The goal of this paper is to 
describe a common dynamic arising in such cases, gen- 
erated solely by kinetic constraints and the geometry of 
the pore and diffusing objects. 

A clear departure from linear diffusion takes plače 
when some object enters a pore, and due to size or geo- 
metrical orientation cannot pass all the way through. If 
there are many such blockers on one side of a membráne, 
carried in a bath of smaller free-flowing particles, we can 
produce a rectifier. Flow in one direction sweeps blockers 
into pores, interrupting further flow. Flow in the other 
direction clears the pores, re-establishing free flow. This 
mechanism was in fact found in the squid giant axon, by 
Armstrong and Binstock in 1965 [2J. 

A tacit assumption of many subscqucnt membráne 
studies is that there must exist forces which hold blockers 
in plače, so that they are not dislodged by thermal fluc- 
tuations. Among proposed mechanisms are electrostatic 
attraction, due to charges fixed in the pore walls [3], 
and microscopic flaps and tethers, to mechanically hold 



a blocker [H [S] . Here we argue that kinetic constraints 
alone can readily perform this function. A blocker must 
backtrack to reopen a pore, and as concentration is in- 
creased, more and more objects obstruct this motion. 
Thus pore clearance can become statistically unlikely, re- 
quiring a serieš of favorable ŕluctuations [5J [7] . 

Figúre 1 illustrates a number of situations where block- 
ing can arise, diffusing objects may arrive at a configura- 
tion that blocks flow through the pore. This is a conse- 
quence of simple geometry only, and can occur in models 
cmploying Hamiltonian or Monte Carlo dynamics, in con- 
tinuous space or on a lattice, with two species of diffusing 
object, or just one. 

The common theme of all of these models is that flow 
through a pore enters a blocked configuration, a dead 
end in the configuration space. In order for the block to 
be removed, objects must back up a finite dištance along 
the path they travelled. If there is a high concentration 
of objects in and around a pore, this retracing may be 
unlikely, and the systém will be trapped in a long-lived 
metastable blocked state. This blocked macroscopic state 
then becomes a type of attractor, in the sense that most 
microscopic configurations lead to it. 

Given a high cnough concentration of particles, a pore 
becomes blocked, and remains blocked for a time that 
depends on concentration and pore geometry, which de- 
termine the strength of the attractor. When the stability 
of the attractor is weak, local density ŕluctuations around 
the pore may cause ŕluctuations between blocked and un- 
blocked pore states, with a consequent intermittent flux 
transmission through the pore. A hard disk simulation 
of a membráne consisting of a single pore with the geom- 
etry oŕ Fig. l(a) displays this behavior [5], as shown in 
Fig. 2. This flux intermittency is qualitatively similar to 
data obtained from patch clamp observations of current 
through a single celí membráne pore [9]. 

The interactions between particles and pores that pro- 
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FIG. 1: Blocking of membráne pores with different pore geometries and different particle types. (a) Large and small disks 
moving in continuous space. (b) Large and small squares moving on a lattice, with a spacing of the size of the small squares. 
(c) Dominos moving on a lattice; while in a pore, a domino cannot rotate. (d) Single species of square, moving on a lattice 
with spacing one-half the size of the square. (e) Single species of disk, in the continuum. In all cases, blockers must move 
upward to re-establish flow. In model simulations, particles in (a) and (e) move in continuous two-dimensional space, governed 
by Hamiltonian dynamics. Particles in (b-d) move on a lattice, with Monte Carlo updates. 



duce the blocking states illustrated in Fig. 1 are simple 
enough that we might expect this situation to occur quite 
generally. We now present an illustrative model that uses 
macroscopic variables to characterize the blocking state 
of the whole membráne, capturing qualitative behavior 
common to all the various blocking scenarios. 

The model membráne contains a large number of pores, 
each of which can either be open or closed. For open 
pores, we assume Fickian diffusion and, for simplicity, 
that one side of the membráne is held at zero concentra- 
tion. The flux across the membráne is thus J = fDC, 
where / is the fraction of open pores, D is a diffusion co- 
efficient, and C the concentration on the occupied side of 
the membráne. We describe the probabilities of a single 
pore entering or leaving the blocked state. 

The probability that during a small time interval a 
blocker enters a pore we take to be proportional to con- 
centration, and a coefficient refiecting the fraction of 
blockers, or blocking configurations, P c ioscs = BC. The 
simulations indicate this assumption to be reasonable. 

For a blocked pore to become open, the blocker must 
experience a fluctuation which backs it out of the pore. 
This will become increasingly unlikely as the concentra- 
tion behind it increases, and as the pore length increases. 
If we assume that a serieš of independent favorable fluc- 
tuations are required for the blocker to move backward a 
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FIG. 2: A "patch clamp" time serieš of the flux of small 
disks through a single pore in the presence of blockers, as in 
the geometry of Fig. l(a). The particles are diffusing through 
the pore from a chamber kept at constant concentration, into 
a vacuum. The flux values are the number of particles passing 
through the pore in each interval of 1000 time steps [8]. 



finite dištance, individual fluctuation probabilities multi- 
ply, and we are led in the simplest case to an exponential, 
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, where C is the concentration difference. 



G is a geometrical factor that encodes the interaction 
between the diffusing blocker and the pore, and would 
depend, for example, on the length of the pore, or the 
dištance the blocker has to retrace. 

Now we consider a chamber of unit volume, leaking 
through a membráne with many pores, into a vacuum, 
and assume the number of pores is large enough so that 
we can treat the fraction of open pores, /, as a continuous 
variable. We obtain a pair of coupled nonlinear first-order 
ODE's. 
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= -fDC, 
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The first describes the change of concentration in the 
chamber due to flux through the available open pores, 
and the second the change in the fraction of open pores 
due to the dynamics as described above. 

If the concentration in the chamber is held constant, 
we can solve for the flow as a function of concentration in 
the steady state, obtaining J(C) = DC/{BCe GC + í). 

Figúre 3 displays the results of a hard-disk effusion sim- 
ulation [10] , along with curves from the above expression 
[TT] , The dotted curves are obtained from a Hamilto- 
nian systém of two sizes of elastic hard disks moving in 
a chamber, one wall of which is a membráne perforated 
with tapered pores with the geometry of Fig. l(a). The 
smaller disks are able to effuse through the pores into a 
vacuum, any that do so are replaced in a random position 
in the chamber, thus keeping concentrations constant. 
We examine four different pore lengths. At a given pore, 
the outward flow may resemble the intermittent behavior 
of Fig. 2, but we take an average over multiple pores, and 
over a long time. 

Plotted is the flux through the membráne as a func- 
tion of the incremental pressure of small disks added to 
the chamber, with a hxed number of large blocker disks 
present. This pressure is the appropriate variable govern- 
ing the dynamical flux in this Hamiltonian systém, and 
is computed from the concentration (disk area) using the 
hard disk equation of state [Ti2] . 
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FIG. 3: Steady-state flow of small disks through a membráne 
pore with the geometry of Fig.[TJa), as a function of pressure 
in the chamber. The dots are data points from a simulation, 
and the curves are generated by the steady-state equation 
above. The relative pore lengths are 2, 3, 5, and 10 (from the 
top curve down). Flows through the reversed pore geometry 
are plotted as negatíve values. 



For low concentrations of small disks, the flow is of 
Fick's law form, with an effective diffusion constant D. 
But as the concentration is increased, large disk blockers 
are swept into pores, and the mean flow is decreased. 
At high concentrations, the probability that a blocker 
can back out of the pore becomes small, and the flow is 
reduced to near zero. 

If the membráne geometry is reversed, that is, disks 
strike the membráne from the Hat side, or bottom-up in 
Fig. l(a), the large disks are not trapped in the pore by 
trailing small disks. In this case, the flux remains a lin- 
ear function of pressure up to high concentrations, these 
values are plotted as negative values in Fig. 3. Note that 
the asymmetric geometry of the membráne lcads directly 
to asymmetric diffusion in the two flow directions. 

The solid curves of Fig. 3 are plotted from the ana- 
lytic expression, we take G values in proportion to the 
simulation pore lengths, and fit the theoretical curves to 
the simulation data with a single overall scaling. The 
blocking factor B varies from 1, when all disks present in 
the chamber are large-disk blockers, to smaller values, as 
an inereasing fraction of small disks are added [TI]. The 
diffusion constant D is found from the simulation data, 
in the linear low concentration range. 

The fit is quite satisfactory, considering the pores in the 
simulation are actually conical, and that density fluctu- 
ations of the hard disk gas are nontrivial [13] . In fact it 
will be more difficult for a blocker to escape than our ex- 
ponential estimate predicts, because a fluctuation which 
empties a pore so that a blocker can back out will ereate 
a locally larger concentration near the mouth of the pore. 
For extended particles, such a fluctuation will be reflected 
in a very large local pressure inerease, and will be highly 
unlikely as concentration inereases [T^l [T3] . The model 



equations will not be particularly aceurate for all pore 
and blocker geometries. But they may be a minimal set, 
deseribing deviations from the linear regime of nonequi- 
librium flow. All of the geometries in Fig. [TJ generate 
qualitatively similar simulation curves. 

This model has the form of a biased random walk, on 
a line with a reflecting barrier on one end (corresponding 
to the blocked end of the pore) and an absorbing bar- 
rier on the other end (the open end of the pore where 
the particle is absorbed by the bulk). The bias is due to 
the mean flow. Previous discussions of diffusion through 
pores have used similar, though unbiased, Brownian walk 
models. However the physical interpretations of the ran- 
dom walk are completely different, for example deseribing 
a blocker tethered outside a pore [TJ], or conformational 
transformations of the pore itself [15]. 

In our simple model, blockers are biased into a pore by 
a concentration difference, that is, a chemical potential. 
We would like to note that a model of identical form is 
produced by considering a single charged blocker moving 
in a pore, biased by an eleetrical potential maintained 
across the membráne. This might be a minimal model 
for a voltage-controlled rectifier, an important funetional 
component of many celí membranes [TJ . Other models of 
voltage control feature charged sensor levers, and other 
mechanical degrees of freedom which must move to re- 
strict ion flow. The model we present here can produce 
voltage control with no additional degrees of freedom, 
just a rigid pore with appropriate geometry. 

Finally we examine the behavior of a leaky chamber 
filled initially at a high concentration, and then allowed 
to decline in time. In Fig. 4 we compare the model 
equations with a simulation of the domino geometry of 
Fig. l(c) [TJ}], though any of the models of Fig. 1 pro- 
duce qualitatively similar results. While the concentra- 
tion gradient across the membráne is large, the leak rate 
is small. But eventually the seal loosens, and there is a 
relatively sudden dumping of the contents of the cham- 
ber. 

For the domino simulations, note the large effect of 
fluctuations. Even in this simulation, with a chamber 
initially containing 35,000 dominos, and a membráne of 
length 4000 with 500 pores, the onset of the high flux 
región is quite variable. In both the simulation and the 
model equations, there is visible a short transient to an 
attracting state where most pores are blocked, from an 
initial condition where all pores are open. 

Clearly, careful studies of the experimental conse- 
quences of this model are required. For example, the 
statistical properties of curves like that of Fig. 2 should 
be compared with real patch clamp data. Another impor- 
tant issue is temperature dependence. The escape prob- 
ability of a blocker held by a kinetic constraint would not 
be expected to have the Arrhenius temperature depen- 
dence of a blocker held in a potential well. 

Our discussion has emphasized the role of kinetic con- 
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FIG. 4: Simulation of the concentration of a leaky chamber 
of dominos, as a function of time. Shown at the top are six 
runs with the same initial concentration, for two pore lengths 
with the geometry of Fig. l(c). Below are curves generated 
by the model coupled equations, for a few values of G. The 
fiux associated with the G=8 curve is plotted as a dotted line. 

straints and geometry in infiuencing flows of extended ob- 
jects, via the simple requirement that a systém back up a 
finite dištance in its configuration space before proceed- 
ing further. This picture of dead-end alleys in configura- 
tion space may be of some generál use in modeling other 
nonequilibium systems. However, geometry alone cannot 
be the whole story. A real membráne pore carries charged 
ions, and is constructed of proteins with charged residues; 
electric fields must play a centrál part. A model includ- 
ing electric charges fixed to pore walls explains observed 
asymmetric diffusion in fabricated conical nanopores |17j . 
Also, living celí membranes may have evolved a number 
of intricate mechanisms for controlling ionic flow. Never- 
theless we suggest that, due to its simplicity, the block- 
ing mechanism described here will often be present, and 
moreover it may have been the dominánt mechanism for 
simpler pores that occurred early in evolutionary history. 
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